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Abstract. We argue that the recently reported Kolmogorov-like magnetic turbulence spectrum in the cool core of the Hydra A 
galaxy cluster can be understood by kinetic energy injection by active galaxies that drives a turbulent non-helical magnetic 
dynamo into its saturated state. Although dramatic differences exist between small-scale dynamo scenarios, their saturated state 
is expected to be similar, as we show for three scenarios: the flux rope dynamo, the fluctuation dynamo, and the explosive 
dynamo. Based on those scenarios, we develop an analytical model of the hydrodynamic and magnetic turbulence in cool cores. 
The model implies magnetic field strengths that fit well with Faraday rotation measurements and minimum energy estimates for 
the sample of cool core clusters having such data available. Predictions for magnetic fields in clusters for which the appropriate 
observational information is still missing, and for yet unobserved quantities like the hydrodynamical turbulence velocity and 
characteristic length-scale are provided. The underlying dynamo models suggest magnetic intermittency and possibly a large- 
scale hydrodynamic viscosity. We conclude that the success of the model to explain the field strength in cool core clusters 
indicates that in general cluster magnetic fields directly reflect hydrodynamical turbulence, also in clusters without cool cores. 
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1. Introduction 



1.1. Magnetic fields in cool cores 



cooling flows - Magnetic Fields - Turbulence - X-rays: galaxies: clusters - inter- 



seems to be able to reproduce the typical field values of clus- 
ters 2 . 



Galaxy clusters contain magnetised plasma on cluster scales, 
' as radio-synchrotron emission of relativistic electrons reveals 
in form of the so-called cluster radio halos. The origin, the 
strength and geometry of these magnetic fields is still a mys- 
tery. 

With th e emer ging of the first theories of magnetic dynam os 
teatcheloJll95(l lKazantsevlll967t Izeldovich et al.lll990l) . it 
was proposed that the cluster magnetic fields are the product 
of turbulence acting on a seed magnetic field. The seed field 
could be a remnant of non-equilibrium processes in the early 
Universe, or due to some weak magnetisation caused by galac- 
tic outflo ws in the form of winds and radio plasma (for a review 
see lWidrowl2002h . 

Initially, wakes of the weakly sonic galaxy motions were 
thought to be the main driver of the intra-cluster turbulence 1 . 
However, it was realised using numerical simulations of large- 
scale structure formation that the violent mergers of galaxy 
clusters are a much more powerful source of turbulence (e.g. 
lTribbldfT9 93). In numerical simulations of magnetic field am- 
plification, the merger-driven turbulence of galaxy clusters 



It is therefore surprising that the strongest magnetic fields 
seem to be located in the centres of clusters which are dy- 
namically the most relaxed, since the last major merger was 
long ago. These cooling flow (now cool core) clusters had 
time to develop a cool, dense central region due to the cool- 
ing instability of optically thin X-ray emission of cooling gas. 
The magnetic fields reported for these cooling flow regions in- 
ferred by Faraday rotation studies were extraordinarily strong 
(up to 50 yuG) compared to the few fiG fields reported for non- 
cooling flow clusters (for rec ent reviews see Caril li & Tavlorl 
l2002HGovoni & Ferettfcool . It was speculated that the strong 
field s could be a resul t of compression in the cooling flow 
dSoker & Sarazinll990l) . 

However, Chandra and XMM observations revealed that 
the standard cooling flow picture, in which unheated gas cools 
down to neutral gas temperatures, must be incorrect, since the 
expected amount of line emission of cold (< 0.3 keV) gas or the 
expected number of stars formed from the condensing gas was 
not detected 3 . Therefore a heat source has to be present which 



1 e.g. IJaffd <1980h: iRoland Jl98lh: iRuzmaikin et alj 
iGoldman & Rephaelil jl99lh : lDe YoungTll992l) 



[ Dolag etai]<1999h:lRoettigeret ai1ll999bl):lbolag et alJ |2002h 
Heckman et al. 1 1989); Fabian et al. 1 1991); Hansen et al. 1 1995); 
lAllenNl995l):IJaffe & BremeJ|l997l):ISmith et alJjl997h:IO'Dea et all 
<1998l): ICrawford et alj Jl999l): iDonahue et all feOOd): lEdgd 1200 ll): 
lOegerle et alJ fcOOll) : Isalome & Combesl J2003I) : lEdge & FraveJ 
<2003l) 
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balances the cooling of the coldest gas. Since the theoretical 
scenario associated with the term cooling flow has been ruled 
out recently, the only observationally motivated term cool core 
is used in the following. 

The energy losses of the cool core have to be balanced by a 
similar energy injection. There have been two main proposals 4 
for the heat source in cool cores: 



(i) Thermal conductivity, which is close to Spitzer's esti- 
mate, and therefore not suppressed by magnetic fields. This 
would allow the inward transport of heat from the hotter 
environmental intra-cluster medium (ICM) 5 . 

(ii) Dissipation of mechanical energy released by the expansion 
and buoyant motion of radio bubbles inflated by the radio 
galaxies which are typically found in the centre of a cool 
core of a galaxy cluster 6 . 



The thermal conduction scenario (i) as the only heat in- 
jection mechanism into cool cores faces severe problems in 
terms of fine-tuning the required energy injection, and explain- 
ing the existence of cold gas clouds which need sufficient insu - 
lation from the keV gas dSokerll2003l iNipofi & Binnevll2004 . 
Too strong conduction would erase the cool core, but too weak 
conduction cannot prevent the cooling catastrophe. Therefore a 
conductively heated cool core should be unstable. 

The scenario (ii) in which the central radio galaxy balances 
the radiative energy losses of the cool core provides fine tuning 
in the form of a self-adapting feedback mechanism: If the radio 
galaxy activity is triggered by cold gas condensing out of the 
cool core onto the central galaxy, the gala xy activity in creases 
until it disrupts further accretion (Chu razov et all200lh . 

In this work we assume scenario (ii), not only since it is - at 
least in our view - theoretically more compelling, but also be- 
cause it predicts a certain level of hydrodynamical turbulence, 
which can be compared to the level required to explain the cool 
core magnetic fields by magnetic dynamo theory. 

There have been some reports on observed signa- 
tures of turbulent flows in cool cores of galaxy clus- 
ters, and the scenario (ii) investigated here seems to be- 



come w 
sis (e.g. 


idely accepted, at least as 


a working hypothe- 


Loewenstein & Fabian 


1990; 


Churazov et allboOll 


Bohringer et all2002UChurazov et alJ2002ll2004. 



4 There are also other mechanisms di scussed in the literature, e. 
the influence of a cosmic ray population I Chandran 2004 ; Cen 200. 
3 Malvshkin 



(2001); Naravan &Medvedev 
J2002I): iRuszkowski & Beeelman 1 2002); 
SokeretalJ <2004: 



Jubelg; 
(2004) 



Voiet & Fabian 
as et alJ l2004h : IChandran & Mar 



(2001) : IVoigtetal 
Cho et alj | 2003t) 
ISokeJ 12004): 



(2004); Soke 
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Dolag etall 



6 Ichurazov et alj bOOll): IBriiggen & Raised f200ll): lOuilis et al 
| 200l|) : IBriiggen et alj <2002h: IChandrart <2004 : iHoeft & Briigger 
(2004); Dennis & Chandran (2005). There have been also observa- 
tions of hot gas bubbles ( Mazzotta et al] |2004 . which may also con 
tribu te to the cool core h eating in a very similar way as the radio bub- 
bles iSoker & Pizzolatol2005l) 



1.2. Observations of cluster magnetic fields 

Magnetic fields in normal clusters and cool core clusters have 
revealed their existence by the diffuse radio halo emission 
in many clusters and radio mini-halos in cool core regions. 
Furthermore, the Faraday rotation of linearly polarised radio 
emission traversing the intra-cluster medium independently 
proves the existence of intra-cluster magnetic fields. If the 
Faraday active medium is external to the source of the po- 
larised emission, one expects the change in polarisation angle 
to be proportional to the squared wavelength. The proportion- 
ality factor is called the rotation measure (RM). This quantity 
can be evaluated in terms of the line of sight integral over the 
product of the electron density and the magnetic field compo- 
nent along the line of sight. 

Magnetic fields in non-cool core clusters of galaxies de- 
tected through Far aday rotat i on me asurements are of the or- 
der of a few fiG. iKim et all dl99lh measured field strengths 
of about 2 fiG on scales of 10 kpc in a statistic al sample of 
point s o urces observed throu gh the Coma cluster. | C larke et all 
(2001), Clarke] J2004 . and Ijohnston-Hollitt & Ekersl J2004 
derive similar field strengths of several fiG in various sam- 
ples of point sources observed within and through various low 
redshift clusters in the northern and southern hemisphere. The 
analysis of RM maps of extende d radio sources have lead to 
the same conclusion (e.g. Feretti et al. 1995, 1999; Taylor et al 



200 ll ICovoni et alJkoOlUF.ilek & Owenll2002UVogt & En Blin 
120031) . 

However, the analysis of Faraday rotation measurements of 
extended radio sources in the centre of cool core clusters reveal 
higher magnetic field strength s. Such an analysis h as been done 
for the Centaurus cluster bv iTavlor et all |2002fc for A 195 8 
(b etter known as 3C29 5) bv lPerlev & TavlorNl99ll). for A1795 
bv lde & Owenl ll 19931). for Cygnus A bvbreher et all i 19871) 
and for Hydra A bv ITavlor & Perlevl Jl993h . From these RM 
measurements, magnetic field strengths of 10 to 40 fiG have 
been reported for the cores of these cool core clusters on scales 
of 3-5 kpc. 

These rather large field values for cool core clusters have 
been revised in the case of Hydra A by a recent analysis 
of the observational data. A high quality Faraday rotation 
map of the north lobe of Hydra A produ ced by the nove l 
PACERMAN algorit hm JDolag et aljEoOst IVogt et al.ll2005l) 
based on the data of ITavlor & Perlevl i 19931) was analysed by 
IVogt & EnBiinl i20()£ . They used a maximum likelihood esti- 
mator for the derivation of the magnetic power spectra, b ased 

on the theory of turbulent Faraday screens (EnBlin & Vogt 

I 1 1 — — I 1 J 

2003; Vogt & EnBlin 2003), and also using the most up-to date 

gas density profile of the cool core, which turned out to make a 

crucial difference. Thereby, a magnetic field strength of l±2fiG 

was found in the centre of the cool core region of the Hydra A 

cluster, which is still a significa ntly large r field than reported 

for non-cool core clusters. Vogt & EnBlin ( 2005) measured the 

detailed magnetic power spectrum from the Hydra A dataset, 

which revealed a Kolmogorov-type spectrum on small scales 

indicating turbulence (see Fig.Q. 

It has been debated whether the magnetic fields seen by the 

Faraday effect exist on cluster scales in the ICM, or in a mix- 
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ing layer around the radio plasma which emits the polarised 
emission jBicknell et all Il990t iRudnick & Blundelll l2003t> . 
However, there is no valid indicat ion of a source loc al Faraday 
effect in the discussed cases (Enfil in et alJ 12*003). and the 
Faraday rotation signal excess of radio sources behind clusters 
compared to a field control sample strongly supp orts the exis- 
tence of strong magnetic fields in the wider ICM (|Clarke et alJ 
l200lUjohnston-HoUitt & Ekersl2004HClarkel2004 ~ 

The morphological structure of magnetic fields in galaxy 
clusters is difficult to obtain due to the projection effects in 
radio observations. However, there are situations in which 
the magnetic fields are illuminated in only sheet-like sub- 
volumes of clusters. This happens whenever short-lived ultra- 
relativistic electrons are injected at a shock wave travelling 
through the ICM. The electrons usually lose their energy be- 
fore the shock wave can travel far away. The electrons thereby 
produce synchrotron radio emission from a nearly sheet-like 
volum e. This emission likely forms the so-called giant radio 
relics fenfilin et alJll998al iRoettiger et alJll999al) . High reso- 
lution radio polarisation maps of the relic in Abell 2256, which 
is one of the largest radio relics known, reveal that the magnetic 
fields are or ganised in filaments or sheets with an aspect ratio 
of at least 5 JClarke & EnBlinl2005l) . Furthermore, the Faraday 
rotati on map of 3C465 reveals stripy patterns ( Eile k & Owenl 
2002]), also suggesting the existence of intermittent fields in 
galaxy clusters. The presence of thermally isolated elongated 
cool H-ff filaments in the core of the Centaurus cluster is also 
best understood by t he existence of filamen tary magnetic fields 
in that environment JCrawford et alJ2005l) . 

1.3. Structure of the paper 

The paper is structured as follows: In Sect.|2]we introduce the 
turbulent dynamo and summarise its expected and observed 
signatures in galaxy clusters, mainly based on the abovemen- 
tioned estimate of the magnetic power-spectrum in the Hydra A 
cluster core. In Sect. [3] we develop a steady-state analytical de- 
scription of the expected magnetic turbulence due to the stirring 
motion of buoyant radio bubbles from the central galaxy, which 
is assumed to regulate the energy content of the cool core. This 
is applied in Sect.[4]to a number of prominent cool core clus- 
ters, and compared to existing information on magnetic fields 
whenever available. Sect.|5]contains our conclusions. 

Throughout the paper, we assume a Hubble constant of 
Hq - 70km/s/Mpc, and translate literature values of lumi- 
nosities, length scales, electron densities and magnetic field 
strengths derived from Faraday rotation measurements to this 
value. 



2. Turbulent magnetic dynamo 

2.1. Dynamo concepts 

The Kolmogorov-like magnetic power spectrum in the cool 
core of the Hydra A cluster indicates that the magnetic fields 
are shape d and probably am plified by hydrodynamical turbu- 
lence (e.g. lDe Young! 1992l) . Therefore, it seems most promis- 
ing to seek the origin of the observed magnetic power spec- 
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Fig. 1. Magnetic turbulence in the centre of the cool core clus- 
ter Hydra A as derived from the Fara day rotation map o f the 
northern radio lobe of Hydra A bv lVoet & Enfilinl ((2005). An 
angle of 45° was assumed in this figure between the line of 
sight and the approaching north-lobe. Variation of this angle 
changes the overall normalisation of the spectrum, but not its 
shape. The right-most data point is likely to be contaminated by 
observational noise in the Faraday map. A central root-mean- 
square magnetic field strength of B lms -13+ 0.2 + 2 yt/G and 
a magnetic aut ocorrelation length of A r = 2.8 + 0.2 + 0.5 kpc 
was derived bv lVogt & Enfilinl ( 120 05). The first errors are the 
statistical uncertainties due to the limited statistics, whereas the 
systematic error reflects the uncertainties in the geometry of the 
radio source and the Faraday screen. 



trum in the theories of turbulent dynamos. A very similar view 
and ap proach are used in the recent work bv lSubramanian et alJ 
(2006), in which the magnetic fields in non-cooling core clus- 
ters were also assumed to be maintained by turbulence. In that 
environment the turbulence is due to merger events and galaxy 
motion, whereas here it is due to the inflation and buoyant mo- 
tion of radio bubbles. 

Since cool cores of galaxy clusters are not believed to ro- 
tate, the gas flow is probably non-helical and the galactic dy- 
namo theories do not apply here. Instead, the non-helical turbu- 
lent dynamo (also called small-scale dynamo) should operate if 
the gas flow is sufficiently random, as it would be in the case of 
developed turbulence. 

The physical details of the small-scale dynamo are still de- 
bated. However, as it turns out, the exact nature of the small 
scale dynamo is of minor importance for an understanding of 
cool core magnetic fields. The observed magnetic field is prob- 
ably determined by the saturated state of such a dynamo, pro- 
vided the dynamo had sufficient time to amplify seed fields to 
the dynamically relevant strength. This condition is fulfilled if 
either the dynamo is very efficient, permitting the very weak 
primordial fields to be amplified, or if the dynamo can start in 
an existing, relatively strong magnetisation, which is only a few 
orders of magnitude below the saturation level. There are argu- 
ments in favour of both pre-conditions being fulfilled, which 
we discuss briefly. 



4 



T. A. EnBlin & C. Vogt: Magnetic turbulence in cool cores 



The random gas motion would stretch and fold any ini- 
tial seed magnetic fields and lead to an exponential growth 
of the magnetic energy density with time with the character- 
istic time-scale being the eddy-turnover time. This proceeds as 
long as the dynamical back-reaction of the magnetic field is 
unimportant. The folding operations of the flow form small- 
scale magneti c reversals perpen dicular to the local field direc- 
tions ( ISchekochihin et alJ2002i) . Magnetic diffusivity limits the 

— 1/2 

scales to be of the order of Ab ~ Aj R m , where Aj is the 
turbulenc e injection scale, and R m is the magnetic Reynolds- 
number (Ruzmaikir^ui2l98s3). The typical bending radius of 
the fields should be of the order of the turbulence length-scale 
At. 

This picture has been criticised by Goldshmidt & Rephaeli 



(1993) as being unable to explain cluster magnetic fields. 
Their main objection is that the transverse size of the mag- 
netic structures Ab should be extremely small, since the mag- 
netic Reynolds number is typically R m ~ 10 28 - 29 in a clus- 
ter environment, leading to Ag of the order of a light second. 
However, magnetic structures with kpc size are required to ac- 
commodate the observed Faraday rotation values with phys- 
ically plausible cluster magne tic fields. This reasoning lead 
iGold shmidt & Rephaeli, Jl993h to argue for a galactic origin of 
the intra-cluster magnetic fields. Rough estimates of the mag- 
netisation of the intra-cluster medium indicate that galactic out- 
flows in magnetised winds and relativistic plasma jets should 
lead to a substantial seed magnetisation 7 . 

The small-scale dynamo picture can only be reconciled 
with the observations if the huge magnetic Reynolds number is 
replaced by a much lower effective value. It was pointed out by 
ISubramanianI i 19991) that the gas motions induced by the mag- 
netic forces lead to a diminishing of the magnetic field strength, 
which can be expressed approximately as an effective mag- 
netic diffusivity. This diffusivity increases with growing mag- 
netic field strength, leading to a decreasing effective magnetic 
Reynolds number permitting the magnetic structures to grow 
to larger spatial scales. On larger scales more turbulent energy 
density is available in a Kolmogorov cascade, thus larger field 
strengths can be accommodated. Larger field strengths imply 
a further decreased effective magnetic Reynolds number and 
therefore a further growth of the fields in strength and length 
scale. This process continues until a saturated state is reached. 
The magnetic e-folding time is of the order of the turbulent 
eddy time scale and might be sufficient to amplify even primor- 
dial magnetic fields in a cluster environment to their observed 
strength. 

It was argued by Sche kochihin et all J2005i) that plasma 
instabilities should lead to an accelerated regime of mag- 
netic amplifications during the pre-saturation phase. Assuming 
an effective description for the plasma particle pitch angle 
scattering by the plasma- waves generated in the instabilities 
Sch ekochihin et all J2005t) showed that a nearly explosive pro- 
duction of magnetic fields from very weak primordial seed 



7 lReed il987i). iDalv & Loebl ll99d). IChakrabarti et alJ <1994l) 
lEnfflin et alJ <1997l Il998bl). iKronberg et alJ jl999l) . IVolk & Atovanl 
( 200Qi). lKronberg et alJl2001l) . lBertone et alJ J2005I) 



fields on cosmologically negligible times-scales might be pos- 
sible. 

Thus, although the details of the generation of dynamically 
relevant magnetic fields in galaxy clusters are still unclear, it 
seems that there are sufficient sources of magnetisation, and 
sufficiently efficient dynamo mechanisms present to amplify 
such fields in cosmologically short times. The question for the 
interpretation of the observational data is in which state do we 
expect the fields to be at present. This requires an examination 
of the saturated dynamo state. 

2.2. Saturated dynamo state 

We start our investigati on of the s aturated state with the view 
introduced by Subra manianl dl999l) . As soon as the Lorentz 
force becomes sufficiently strong, the fields do not follow the 
flow passively, but try to disentangle themselves. This back- 
reaction motion leads to an increase in the effective magnetic 
diffusivity, and therefore to a lower, renormalised magnetic 
Reynolds number implying that fields become organised on 
larger scales. The Reynolds number decreases until it reaches a 
critical value R c ~ 10 - 100, below which the turbulent dynamo 
would stop operating. The system reaches a saturated state, in 
which the magnetic correlation length scale 



At R. 



-1/2 



(1) 



is solely determined by the turbulence length scale and the crit- 
ical _Re^noJ^sjnimber. 

ISchekochih in & Cowlevl £2006 ) acknowledge that their ex- 
plosive dynamo alone fails to explain the observed level of 
cluster magnetisation, but only by up to one order of magni- 
tude. Therefore, the remaining amplification and the saturated 
state should result from a conventio nal non-helical dyn amo 
(e.g. ISubramanianI 1 1 9991) . However. ISchekochihin & Cowlevl 
(2006) speculate that the physics of the plasma instabil- 
ity driven dynamo still imprints on the saturated state, 
leading to a characteristic field strength with Ab ~ 

(VthjiA'T) (pi/ ^t) • Here, v t h and vj are the ion and the 
turbulent velocity, respectively, p; is the gyroradius of an ion 
within a magnetic field that is in equipartition with the turbu- 
lent energy density. This picture of the saturat ed state can be 
mathematically mapped onto the model of Subramanian (1999) 
if we define an effective magnetic Reynolds number of 



R 



fth,i 



(2) 



which takes values in the range 10 2 3 for cl uster environments . 

In the pictur e initially developed b y iBatcheloJ {l950); 
Kazantsev ( ll967l) : Zeldovich et al l dl990l) . the magnetic fields 
are highly intermittent: only a small fraction f B of the volume 
is actually strongly magnetised, whereas the remaining vol- 
ume 1 — fn does not carry dy namically imp ortant magnetic 
fields . Iz eldovic h et alJ Jl990|) and others (|Ruzmaikin et alJ 
ll989HSokolovetalJll990l: ISubramanianI 1 19991) assumed the 
magnetised regions to be organised as magnetic flux ropes 
of volume At A 2 b per eddy volume A T . The occupied volume 
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fraction would be fg ~ A 2 B /A T ~ Z?" 1 . Such highly filamen- 
tary structures were not confirm ed by numerical simulations. 
By using a spectral MHD code, ICho & Vishniad (|2p00) find 
for example that the magnetic autocorrelation functions exhibit 
only a moderate anisotropy on scales where the magnetic en- 
ergy density peaks. However, numerical simulations in physi- 
cal space indicate that the magne tic structures are more sheet- 
like with thickness Ar ferandenburg & S ubramanian 2004; 
ISchekochihin & Cowlevll2006h. leading to a f g ~ An I At ~ 
R^ 1 ^ 2 as assumed in Subramani an et al] J2006h . Furthermore, 
as we explained in Sect. 11.21 radio observations of the cluster 
radio relic in Abell 2256 also strongly support the assumption 
of intermittent magnetic fields, probably of sheet-like structure. 

In order to allow for any geometry of the magnetic struc- 
tures in the saturated state, we assume that they are effectively 
c/-dimensional, giving 



fs 



A T A B 



Rr 



with 1 < d < 3. 



(3) 



Under the steady-state conditions of the saturated dynamo, 
the hydrodynamical dissipation of turbulent energy on scale 
At and the relaxation of magnetic structures bent on the same 
scales Aj should have the same time scales: 



tb 



t t 



Vt 



9 „struct. 



P 

2 ex ' 



(4) 
(5) 



where v^ ract 



and sf uct - 



are the Alfvenic velocity and energy 
density within the magnetic structures, and Vj and sj the tur- 
bulent velocity and energy density, respectively. Then Tg ~ tj 
implies that the magnetic fields within the structures are in 
equipartition with the environmental turbulent energy density, 
and the volume-averaged magnetic energy density is therefore 
lower than these by the magnetic volume filling factor: 



f B ~ s T fg ~ e T R c 



(6) 



Thus, with the knowledge or assumption of the critical 
Reynolds number R c , and the effective dimensionality d of the 
magnetised regions it is possible to translate properties of the 
hydrodynamical turbulence to the magnetic turbulence and vice 
versa, under the assumption that the system is in the saturated 
dynamo state. 

In the numerical examples, we will investigate three sce- 
narios, which differ in assumed magnetic topology and critical 
magnetic Reynolds number: 

1. Magnetic fields mostly organise in flux tubes (d = 1, 
R c = 20, => fg = 0.05). This scenario reflects the 
original assumption o n turbulent magnetic structures by 
Zeldovich et al but modified with the concept o f 

the renormalised Reynolds number ( Subr amanianl ll999). 
Here we consider it to investigate the induced hydrody- 
namical viscosity on large scales, since there are pub- 
lished predicti ons for t his quantity in the flux rope scenario 
( Longcope_£Lal|20J}3j). 



2. Magnetic fields mostly organise in sheets and ribbons (d = 
2, R c = 35, => f B = 0.17). This scenario is more in 
agreemen t with recent nu merical simulations, as argued by 
ISubramanian et al.ld2006l) . 8 

3. The assumed saturated state o f the explosive dynamo of 
ISchekochihin & Cowlevl (|2006), which we characterise by 
d — 2 dimensions and by an effective Reynolds number 
R c - /?* provided by Eq.|2] 

The scaling of any quantity with the parameters R c and d will 
also be provided, so that this theory may be applicable even if 
the adopted values need readjustment. 

2.3. Magnetic viscosity 

In the saturated dynamo state, the magnetic fields do not pas- 
sively follow the hydrodynamical flow of the bulk motion, but 
possess an independent velocity component due to the Lorentz 
force. This should lead to a slippage of strongly and weakly 
magnetised regions relative to each other, possibly causing 
some friction due to induction of small scale eddies in the 
wakes o f the magnetic structur es. For the flux rope dynamo 
scenario, Longcop e et al 1 ll2003l) give the expected viscosity on 
large spatial scales in the saturated state to be of the order of 
4% of the turbulent diffusivity kt ~ Vt /It/3, thus 



0.04 



(7) 



Also in the scenario of mainly 2-dimensional magnetic struc- 
tures, one would expect such viscous effects, but no estimates 
of the viscosity exist to our knowledge. 

2.4. Signatures 

According to the scenario described above, the magnetic fields 
in the turbulent cool cores of galaxy clusters should exhibit the 
following properties: 

A. The magnetic power spectrum should either reflect the hy- 
drodynamical power spectrum due to their dynamical cou- 
pling, just with a lower normalisation and on smaller spa- 
tial scales or the small-scale magnetic fluctuations are ex- 
pected to follow a Goldreich-Sridar law, which is indis- 
tinguishab le from a Kolmogorov spectrum in the isotropic 
average (Gol dreich & Sridhar1ll997l) . That means that a 



8 lHaugen et alJ 120041) numerically investigate magnetic turbulence 
in the case of a Prandtl number of order unity, and find a magnetic to 
turbulent energy ratio of the order of f B ss 0.4 in this regime. Since 
they report a R c a 35, their case could be describe d by d = 2.5 di- 
mensi onal magnetic intermittency. However, as Dennis & Chandran 
12005) point out, the dissipation rate of turbulent energy seems to be 
doubled compared to not (or less) magnetised scenarios, which would 
lead to a reduction in the turbulent and magnetic energy densities by 
a similar factor. For the observables of the magnetic fields, like field 
strength, length and Faraday signal, there would be very little numer- 
ical difference to the d = 2 scenario. For this reason, we do not follow 
this scenario separately, but assume it to be subsumed under the d = 2 
case. 
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Kolmogorov-inertial range power law behaviour of hydro- 
dynamical turbulence should also be found in the magnetic 
spectrum in any scenario. 

B. The average magnetic energy density s B is lower than the 
turbulent kinetic energy density e^n by Sb ~ £kin /b, where 
/b ~ R^ d ^ 2 ~ 0.05 . . . 0.2 depending on the critical mag- 
netic Reynolds number and the dimensionality of the mag- 
netic structures. 

C. The magnetic fluctuations are concentrated on a perpendic- 
ular scale Ab, which is smaller than the hydrodynamical tur- 
bulence injection scale Aj by Ab ~ AjR^ 1 ^ 2 . 

D. Magnetic correlations exist up to a scale Aj, turn there into 
an anti-correlation as a consequence of V ■ B = 0, and 
quickly decay on larger scales. 

E. The fields may be spatially intermittent. This may be under- 
stood by Zeldovich's flux rope model, in which magnetic 
ropes with diameter Ab are bent on scales of the order Aj. 
Alternatively, as also supported by recent numerical simu- 
lations, it can be understood in terms of magnetic sheets of 
thickness As and size /It- 

F. In the case of intermittent fields, the field strength within 
the magnetic structures should be in energy equipartition 
with the average turbulent kinetic energy density of their 
environment. 

G. The magnetic drag of such intermittent structures produces 
a hydrodynamical viscosity on large scales. In case of the 
Zeldovich-flux rope scenario, an estimate of a viscosity of 
4% of the turbu lent diffusivity kj ~ Vx Aj/3 was made by 
iLonecope et alJd2003l) . 

2.5. Observations 

We investigate briefly if the above predictions of the non- 
helical dynamo theory are in agreement with observations. 

A. A Kolmogorov-like magnetic power spectrum in a clus- 
ter cool core is revealed by the Faraday rotation map of 
Hydra A JVogt & Enfilinl2005l) . 

B. Translating the Faraday rotation based estimate of the 
magnetic energy density to the expected turbulent en- 
ergy density sj ~ sb//b in the Hydra A cluster yields 
0.3... 1.0 • 10~ 1() erg cirT 3 , which corresponds to turbu- 
lent velocities of v tur b ~ 250 . . .430km/s. This is com- 
parable to velocities of buoyant radio plasma bubbles 
(EnBlin & Heinz 2002), which are e xpected to stir up tur- 
bulence (e. s. IChurazov et all200ll) . 

C. The expected turbulence injection scale in the Hydra A 

1 /2 

cluster core is of the order of A T ~ A B R^' ~ 10 ... 20 kpc, 
again consistent with the radio plasma of Hydra A being the 
source of turbulence, since the turbulence injection scale 
and the radio lobes of Hydra A have comparable dimen- 
sions. The dynamical connection of the radio source length 
scale and the magnetic turbulence scale would explain why 
the Faraday map of Hydra A is conveniently sized to show 
us the peak of the magnetic power spectrum (see Fig.0. 

D. The expectation of weak magnetic fluctuations on scales 
larger than the turbulence injection scale is hard to test with 
the available data due to the limited size of the RM map 



used. However, the downturn of the magnetic power spec- 
trum at small fc-values visible in Fig.^is in good agreement 
with the requirement of magnetic anti-correlations on large 
scales. 

E. Magnetic structures in the form of flux ropes or ribbons 
might have be en detected as stripe d patterns in the RM 
map of 3C465 feilek & Owenl2002i) . and as polarised syn- 
chrotron filaments in th e cluster radio relic in Abell 2256 
llClarkefeFnminl E005). in support of the assumed inter- 
mittence of the small-scale dynamo . 

F. The fraction of the strongly magnetised volume may be 
as small as /b = R^ 1 ~ 0.05 . . .0.2. Strongly intermittent 
fields in galaxy clusters may help to reconcile the discrep- 
ancy between Fara day and inverse Co mpton based mag- 
netic field estimates (EnBlin et al. Il999i) . The latter estimate 
could easily be biased to low field regions due to the faster 
removal of ultra-relativistic electrons in strong field regions 
by synchrotron emission. 

G. The expected hydrodynamical viscosity on large scales in 
the Hydra cluster is of the order of 10 28 cm 2 /s (if the esti- 
mate in Eq .|7]based on the flux rope picture is applicable). 
Fabi an et alJ ( l2003bl) argue, for the comparable Perseus 
cluster cool core, for a lower limit on the large-scale vis- 
cosity of 4 • 10 27 cm 2 /s. They base their arguments on the 
observation of very elongated H-a filaments which suggest 
a laminar flow pattern behind buoyantly rising radio bub- 
bles. An upper limit on the viscosity in the (non-cooling 
core) Coma cluster o f ~ 3 ■ 10 29 cm 2 /s is reported by 
ISchuecker et alJ J2004I) . Both limits are consistent with our 
coarse estimate of the large scale viscosity and enclose it. 
However, the presence of a significant viscosity, as well as 
the underlying magnetic intermittency, in galaxy clusters 
should be regarded as being speculative, and not directly 
confirmed by observations. Nevertheless, there is a grow- 
ing number of theoretical investigations 9 on the possibility 
that viscosity helps to explain properties of the cluster gas. 

3. Magneto-hydrodynamic turbulence in cool 
cores 

3.1. Scope of the approach 

A theoretical model of the magneto hydrodynamic turbulence 
in cool cores of galaxy clusters is developed in the following. 
The model is intended to capture the most essential features 
of the physical picture. We do not attempt to make accurate 
numerical predictions, but hope to get insight into the scaling 
relations of the different quantities and their order of magnitude 
values. 

Although cool cores exhibit a structured gas distribution, 
for our order-of-magnitude calculation we describe them as 
quasi-homogeneous spheres. We briefly present a radial exten- 
sion of our model in Sect. l3.7l 
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The episodic injection of radio plasma by the central galaxy 
into cool cores should lead to temporal variations of the state of 
the core. Nevertheless, we treat the system as being in a steady 
state. Our estimates provide therefore only approximate tempo- 
ral mean values. Observations of cool cores will always show 
snapshots of the cool core evolution. Thus, deviations between 
our predictions and the observed state of individual cool cores 
should not be too surprising. 

Detailed numerical modelling of the processes would be re- 
quired to overcome our simplifying assumptions which is well 
beyond the scope of this initial investigation of the scenario 
presented. 

3.2. A simplified cool core description 

A cool core is a condensation of cold gas of mass M cc , which 
dropped out of the hot phase of a galaxy cluster due to the faster 
radiative cooling of denser gas. The observables of the cool 
core, which can be used as diagnostics of the physical param- 
eters, are the bolometric X-ray luminosity of the cool core L cc , 
its temperature T cc , and its radius r cc from which the central 
electron density « cc can be deduced. 

The energy feedback of the central radio source prevents 
the cool core gas from falling below a characteristic temper- 
ature 10 T cc ~keV, which is given by the requirement that the 
atomic line emission does not dominate. Otherwise, in the case 
of a lower temperature, the gas in the cool core would rapidly 
condense onto the central galaxy due to a cooling catastrophe 
driven by emission-lines. The pressure in a cool core is 



^cc — 2 ^cc ^cc ? 



(8) 



where we ignore here and in the following the presence of any 
element heavier than hydrogen in pressure and mass terms for 
simplicity. The cool core gas mass is therefore 



An , 

M cc = n cc ntp V cc , where V cc = — r cc . 



(9) 



The X-ray luminosity of the cool core is given by 

^cc = A x (r cc )n^Vc C , (10) 

where 



Ax(r cc ) = A Tl + A Mnes (r cc ) 



(11) 



consists of a Bremsstrahlung- and a line emission term. 
The line emission term should be at best comparable to the 
Bremsstrahlung-term, since otherwise the full core would run 
into a cooling catastrophe. Therefore, we can ignore the line- 
term in our rough estimate 



A X (T CC ) ~ A 

-^CC » 



(12) 



with Ao = 5.96 • 10 24 ergs 1 cm 3 keV (assuming a metal- 
icity of 0.3 solar). 



10 We express temperatures in terms of energies by setting the 
Boltzman constant to unity (fc B = 1) 



3.3. Hydrodynamical turbulence injection 

We assume that the injection from radio galaxies is the domi- 
nant heating mechanism. Although radio galaxies inject energy, 
not all of it is in the form of turbulence. During the initial phase 
of inflation of a radio bubble there are shocks, which c an heat 
the environment ( Heinz et alll998tlFabian et all2003al) . Later, 
a series of sound waves are produced, which might be dissi- 
pated through viscosity in the medium 11 . When a radio bubble 
raises buoyantly in the cool core atmosphere, the environmen- 
tal gas flows around it, which leads to injection of kinetic en- 
ergy from the central radio galaxy into the cool core. The time 
averaged turbulence power of the source is 



Lr - /t L cc - j]j L n 



(13) 



where fj is the fraction of cool core heating by turbulence, 
tjt < 1 is the efficiency of kinetic energy transfer between the 
injected radio plasma and the cool core gas, and L rg is the me- 
chanical lumino sity o f the radio galaxy (PdV work per time). 
IChurazov et alJ(E00 2) argue that the efficiency factor rji should 
be of the order of one since the radio plasma loses most of its 
energy by adiabatic expansion during its buoyant rise through 
one scale-height of the cluster atmosphere. This energy is trans- 
fered mechanically to the kinetic energy of the ICM gas and 
is finally dissipated as heat. It is controversial whether most 
of this energy is dissipated within the cool core, or within the 
outer regions of the cluster. Here, we assume the former and 
set rjj = fj = 0.5 also allowing non-turbulent heating (shock 
waves, sound waves, heat transport). However, we also provide 
the scaling of our results with these parameters. 

The injected kinetic energy is dissipated through a 
Kolmogorov cascade within an eddy turnover-time 



Tj ~ Af/\j , 



(14) 



where Aj and Vj are the turbulence injection scale and turbu- 
lence root mean square velocity, respectively. The average tur- 
bulent energy of the cool core Ej = M cc v 2 /2 is therefore 



Ej - Lj Tj as /t L cc Aj/\j , 



yielding 

2 /t L c , 



'T 

Aj 



2 /t L ct 



,Vc, 



(15) 



(16) 



This would allow the determination of the turbulent velocity 
and energy density if the turbulent injection scale were known. 
In the following, we attempt a rough estimate of the expected 
turbulent length scale Aj. 

3.4. Hydrodynamical turbulence injection scale 

In the picture adopted for this work, the turbulence is stirred by 
the movement of buoyant radio plasma bubbles. The turbulence 
injection scale should therefore be of the order of the radius rbub 



IPringlej Jl989h : IChurazov et alj <200ll): iReynolds et all 120021): 
iFuiita et alJ h.004 ): iRuszkowski et alj J2004) : iDennis & Chandranl 
<2005l) 
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of the bubbles, which we approximate to be spheres. Assuming 
an ultra-relativistic equation of state for the radio plasma, we 
find the volume of the bubble as 



Vbub = £ b ub/(4P), 
where 



(17) 



£bub = L r g Tbub /2 = /t^cc Tbub / (2 77t) ( 1 8) 

is the mechanical energy released by the radio galaxy into a 
bubble during the time Tbub the bubble needs to leave the cool 
core buoyantly. If the turbulence stirring bubbles were filled 
mostly by thermal gas, the volume of the bubble would be 
Vbub = 2 £bub/(5 P), which would make some difference to our 
estimates, but none that change the order of magnitude of our 
results. 

The time available for the jets to inflate the bubble is its rise 
time through the cool core 



T bub 



fee 1 2 F cc 

r bu b/ 2c s ' 



(19) 



and can be estima ted from the ba lance of drag and buoyance 
forces (e.g. EnBlin & Heinz 2002). 

Combining Eq.^lED and^]yields the bubble radius 



( 3/tZ-o 



J/2 \ 



r b ub= ,/ p , (20) 

^647T77TC S f cc ) 

and thereby the required hydrodynamical turbulence scale Aj ~ 
rbub- Using this and the derived relations above, the turbulent 
velocity is 



Vt 



i 2 1 1 _L 
^0 ^ cc ^" cc ^ cc 



2a 52i / T 3 ml 
173km/sl— j ( aicm _ 3 ) ( 10kpc ) (^") ' 



(21) 



J T ) VO.lcm- 3 / \10kpc 
and thus, the turbulent energy density is given by 



3 13 6 



£ T = 



3 21 /^j, tVq tflp Tl cc T cc T cc 

2S5A/J 



Tabled contains the values for Vt, Aj = r^b, and sj as for 
a number of prominent cool core clusters, expected from our 
steady-state description of cool core turbulence. Deviations due 
to episodic evolution of cool core turbulence are possible and 
expected. 

3.5. Magnetic turbulence 

The energy density of the magnetic turbulence eg is lower by a 
factor fs ~ 0.05 . . . 0.2 than the kinematic one ex- The length- 

— 1/2 ^ 

scale is also smaller by a factor R c . The root-mean-square 
magnetic field strength in the cool core is therefore given by 



fin 



= ^7TS T f B = 



2* 3**2 771 Ajm" n£rlT£ 



(23) 



52T/ T 3 tfc 4 

(5.6^G ft)* G^)*(,fc)' 



in our flux rope (d — 1) and magnetic sheet (d = 2) scenarios, 
respectively. The magnetic autocorrelation length is 



Ab - — - 



3 7 Ag ml "Jc^c 



(24) 



= - 53kpC (oTc^) 7 (l^) (keV) 7 ^) 



3.6. Predicting an observable: RM dispersion 

Knowing the depth of the Faraday screen r cc , the magnetic field 
strength (B 2 ) and the magnetic autocorrelation length Ab allows 
one to predict the expected dispersion of the Faraday rotation 
measure from a central radio source 12 , assuming an isotropic 
di stribution of ma gnetic field strengths with the help of Eq. 40 
of lEnBlin & Vogtl(l2003h : 



(RM 2 ) 



1 2 2 

- An IT 



-*() "cc 



r C cA B (B 2 ) 



(25) 



1 n fl^ r/j Aq ml 



(26) 



2^1 



T 1 



rmTZ 



5-/, A', 

f 739 (m\\ I « cc \n 1 r cc \t I t x \—a (R^Y? 

I m 2 \f T ) lo.lcm-'j UOkpc; \keV ) \2Q J 
1 1 29 11 _ 1 _ 1 

1180 />t\3 / \u / rg. \T(Tg,\ W (RA ~- 

V m 2 \f T ) ^O.lcm-V UOkpcj \keV ) \35 ) 



The two cases correspond to filamentary and sheet-like mag- 
netic structure scenarios, respectively. Here, «o = / (2n m 2 c 4 ) 
is the usual Faraday rotation constant 13 . 

The scaling of the Faraday dispersion, which can be writ- 
ten approximately as RM^ S oc n 2 c r^J 2 T^ c , is sufficiently simi- 
lar to that of the bolometric luminosity of the cool core L cc cc 
and also to the mass deposition rate M oc L cc /T cc oc 



1 r i T 

cc ' cc 1 cc 



1/2 
c 

-1/2, 



, Vc , , to expect significant correlations, which are indeed 
observed (e.g. Tayl or et al .120021) . 



(22) 3.7. Approximate treatment of the cool core structure 



In the following, we generalise the simplified geometry of a 
homogeneous cool core to one which has a radial structure. 
This generalisation allows us to predict the radial dependence 



12 Background sources should have twice the variance (RM 2 ) given 
here due to the doubled screen depth. 

13 Eq. l25l differs from the usually used, but inaccurate RM dispersion 
formula, which is based on the cell model for the magnetic fields 



(RM 1 ) * - o£ ;/ 



"cc 'cc arm 



(B 1 ) 



The proper magnetic autocorrelation length A B is used in Eq. I25I which 
differs in a non-trivial way - since this is magnetic power-spectrum de- 
pendent - from the RM autocorrelation length Arm derived from RM 
maps. Usually, A RM > A B . Second, the numerical factor | in Eg. 1251 
properly takes account of the effect of the constraint VB = 0, whereas 
for the above Eq., uncorrelated patches with internally constant mag- 
netic fields were assumed, which do not have V ■ B = at the patch 
boundaries. Due to the differences, published magnetic field values 
based on the above Eq. do not need be accurate, although the two er- 
rors partly compensate each other. 
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Fig. 2. Radial profile of the expected RM dispersion for a line 
of sight starting in the midplane of a fi-profile cool core, but 
displaced from the cluster centre by x ± . The RM dispersion is 
given in units of our result RM^ S in Eq.[25]for the solid sphere 
model of the cool core. 

of the various quantities and to test if our predictions of the 
central RM dispersion are expected to change significantly due 
to the contribution from radii outside the cool core radius. 

In order to provide a radius dependent geometry, we make 
an additional assumption about the radial balance between 
heating and cooling. We assume that the cool core is settled 
in an hydrodynamical state, where heating and cooling are bal- 
anced locally at all radii below the cooling radius. Note that 
this does not necessarily need to be true, since in a convecting 
system the location of heat deposition and radiative cooling can 
differ, but for the sake of simplicity, we use this approximation. 

The cool core is divided into spherical shells of thickness 
dr, which have a volume dV(r) = An r 2 dr and a luminosity 
dL(f) = Ax(T cc ) n 2 ,. dV(r). Our arguments about the balance of 
turbulent heating and radiative cooling made in Sect. l3.3l can be 
identically applied to each of these shells individually, leading 
to a relation equivalent to Eq .1161 



2/ t dL 



dV 



(r). 



(27) 



To derive the local turbulence injection scale, we use again 
the local characteristic buoyancy time T^, u \,(r) ~ r/v^ir) to 
predict the bubble size at any cluster radius r > r cc . For r < r cc , 
we use r cc to treat the central flat density profile properly. In 
doing so, we find that the replacement 



max(r cc , r) 



in any former equation describing local properties as in 
Eqs. 1201231 gives the appropriate function of radius r. Of spe- 
cial interest may be the expected radial scaling of the typical 
magnetic field strength: fi ms oc n]?^\r) r 3 ^ 7 which is usually 
a decreasing function of the radius due to the steepness of the 
electron profiles. For example for « e oc r~ 2 (and T(r) « const) 
we find B lms oc r~ 10/7 oc r^P which lies well within the range 



of usually assumed magnetic field scaling with electron density 
firms oc n° 5 - 1 (e.g. Dola g et all200ll) . 

Global quantities like the energy content, the luminosity of 
a volume, etc. are calculated by performing the appropriate vol- 
ume integrals. The increase of the Faraday dispersion along the 
line of sight should therefore be given by the differential anal- 
ogy to Ea.l25l 



d{RM ) 
~dl 



-a l n l cc A B (B ) 

5 2 - - — 15 

3 21 7r «q 77^ Ag OTp n e 7 (r) max T (r cc ,r) 

70 5 2 2-- - 



(29) 
(30) 



which reduces to the former result given by Eq.|5S]if integrated 
over a cool core sphere with constant properties within r cc . 

Eq.|30|allows us to check under which conditions our con- 
stant core model gives appropriate results, and under which 
conditions the RM contributions from larger radii become es- 
sential. 

For an electron density profile described by the usual 
model, 



n e (r) 



r 2 

1 + - 



(31) 



we find that the contribution to the RM dispersion per loga- 
rithmic radius peaks around the core radius r cc if 5 > 0.25 (a 
similar statement can be made for the cool core dominating 
the X-ray emissivity if /3 > 0.5). Since typical cluster electron 
density profiles are described by /3 > 0.5, our ignorance of con- 
tributions of outer regions to the RM dispersion only introduces 
a moderate error of the order of 20 - 40% (depending on p) for 
lines of sight starting in the cluster centre. Off-centre lines of 
sight, which often occur for extended radio lobes in our sam- 
ple, will result in significantly lower RM dispersions. This can 
be seen in Fig.|2l where numerically estimated /?M-profiles are 
shown. 

Typical pairs of radio lobes are not located in the midplane 
of the cluster. One lobe has a shorter and one has a longer line- 
of-sight through the Faraday-active medium. In a spherical ge- 
ometry of the cluster gas, the additionally accumulated {RM 2 ) 
of the back lobe, and the missing of the front lobe should com- 
pensate roughly in an average over both lobes, provided both 
lobes are in mirror symmetric positions with respect to the clus- 
ter centre. Therefore, in such a case the midplane RM disper- 
sion is a good reference point. 



(28) 4. Application to cool cores 



In the following, we apply our model to a number of cool core 
clusters. The input parameters of our calculations are the cen- 
tral electron density n cc , the central gas temperature T cc , and 
the core radius r cc of the cool core. The numbers are taken from 
the literature and the corresponding references are given in the 
sections refering to individual clusters in the Appendix lAl to- 
gether with detailed discussions of the individual datasets. The 
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Fig. 3. Comparison of observed (cross) and theoretically ex- 
pected (case d — 1, R c = 20) (dashed line) root-mean-squared 
RM signal. The expected RM is calculated for a radio source 
located at the centre of the cluster, whereas most RM measure- 
ments are somewhat off-centre, probably leading to a reduced 
observed signal. 

cluster parameters and the derived properties are summarised 
in TableG] 

As a consistency check for the radii of the cool cores used, 
we calculate the emissivity L cc within the cool core accord- 
ing to Eq. EI| an d compare this to the reported bolometric X- 
ray luminosity of the complete cluster L c i ust er. The fraction 
of X-ray luminosity due to the cool core is in the range of 
^cc/Lciuster = 1% - 50%. The typical fractional luminosity 
within the cooling radius, which is defined by the gas having 
a cooling t i me le ss than the Hubble time, was estimated by 
IPeres et alJ |l?98). Since the cooling radius is larger than the 
core radius of the cool core, a systematic difference is expected, 
and indee d found sinc e L(< r ^piing) /L c i uster = 40% - 60% ac- 
cording to IPeres et alJ dl998l) . The large differences for some 
clusters between the ratios are due to shallow electron density 
profiles with /3 ~ 0.5 (see Sect.E). 

The other quantities in Table ^ are calculated according to 
the formulae given in Sect. [5] Note that our expected RM 1 ^ 
is calculated for a polarised synchrotron source in the middle 
of the cool core, whereas the real radio emitting volume may 
be displaced due to an inclination between radio jet and line- 
of-sight, and/or due to a non-central position. This will cause 
some deviation of our expectations from the observations and 
usually biases the observational values to be lower. 

For a sample of cool cores with luminosities between 
10 43 - 10 45 erg/s, we predict turbulent velocities in the range 
100- 300km/s, magnetic field strengths in the range 3-13 fiG 
(1-d scenario) or 6-23 yuG (2-d scenario), and Faraday disper- 
sions in the range 100 - 4000 irT 2 (150 - 6000 in the 2 - d sce- 
nario). These values should be compared to existing Faraday 
rotation measurements, field strength estimates, and future X- 
ray-spectroscopically determined velocity dispersions, as done 
in Tabled an d Figs.[3]&E] We further list our expectations 
for the large-scale magnetic viscosity, which is speculative, and 



Fig. 4. Expected (case d — 1, R c = 20) and observed RM dis- 
persion values versus the bolometric cluster luminosity. A lin- 
ear correlation is also shown to guide the eye. 

which were estimated using the calculations of lLongcope et alJ 
d2003l) for an assumed 1 - d flux-rope picture of the magnetic 
field configuration. 

To summarise our results, we find that both the one and two 
dimensional scenarios for the magnetic structures predict rota- 
tion measures which are of the same order of magnitude as the 
observed ones. In cases where observational estimates for field 
strength and length scales existed, there is also a better than 
order of magnitude agreement. The explosive dynamo seems 
to underpredict the observed Faraday rotation by a factor of 
two, due to the shorter magnetic length scale predicted. In gen- 
eral, we would not expect our model to be more accurate than 
to within a factor of two. This is especially true for our rough 
parametrisation of the saturation state of the explosive dynamo. 

5. Conclusion 

We showed that many properties of magnetic field measure- 
ments in cool core clusters, and especially in the case of 
Hydra A, which we investigate in detail, seem to support the 
picture that a saturated small-scale turbulent dynamo is main- 
taining the magnetic fields. Although there exist dramatic dif- 
ferences between small-scale dynamo scenarios, their saturated 
state might be quite similar. If it can be assumed to be mostly a 
balance between the turbulent entangling and Alfvenic straight- 
ening of magnetic field lines, the state is characterised by 
only two numbers: the effective dimensionality of the mag- 
netic structures and the ratio between turbulence injection scale 
and magnetic correlation length. Casting three representative 
small-scale dynamo descriptions into this formulation, we find 
that they are all in rough agreement with the data, which is 
as much as we can expect given the approximate nature of 
our treatm ent. The investigated scenarios are the flux rope dy- 
namo (e.g Zeldovich et al. 1990; Subramanian 1999, with a 
renormalised Reynolds numb er), the fluctuation dynamo as de- 
scribed in Subramanian et al. (2006)," and the explosive dynamo 
introduced by Schekochihin et al. (2005). 
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11.8 


8.16 


3.28 


0.14 


L*[*^- '"cooling )/ ^'cluster 


0.52 


0.41 


0.49 




0.61 




0.59 


0.3 


0.47 


0.5 


Lj cc j ^cluster 


0.37 


0.14 


0.3 


0.2 


0.14 


0.0045 


0.51 


0.18 


0.078 


0.005 


hydrodynamical turbulence 




Ab ~ ''bub [kpc] 


7.68 


1.45 


1.7 


2.84 


7.52 


0.2 


13.3 


7.15 


5.94 


0.25 


ej [10~ 10 erg/ccm] 


0.29 


0.16 


0.76 


1.27 


0.35 


0.12 


0.37 


0.13 


0.089 


0.11 


Vt [km/s] 


250 


156 


244 


283 


240 


87.5 


295 


225 


177 


94 


magnetic turbulence 


d = 1 R c 


= 20 


-flux rope dynamo 














As [kpc] 


1.72 


0.32 


0.38 


0.64 


1.68 


0.045 


2.97 


1.6 


1.33 


0.056 


firms b"G] 


6.06 


4.55 


9.78 


12.6 


6.68 


3.8 


6.82 


4.04 


3.34 


3.74 


*v, S c [10 27 cm 2 /s] 


7.89 


0.93 


1.71 


3.31 


7.43 


0.073 


16.1 


6.61 


4.33 


0.097 


RM™s [nT 2 ] 


1520 


351 


1730 


4010 


1930 


91.5 


2590 


607 


400 


96.9 


RM^ s /R]VCs 


0.89 


1.88 


0.69 


0.72 


0.56 


0.27 










magnetic turbulence 


d = 2 R c 


= 35 


-fluctuation dynamo 














A B [kpc] 


1.3 


0.25 


0.29 


0.48 


1.27 


0.034 


2.25 


1.21 


1 


0.042 


firms UlG] 


11.1 


8.36 


18 


23.2 


12.3 


6.99 


12.5 


7.44 


6.13 


6.88 


RM™s [nT 2 ] 


2430 


561 


2770 


6420 


3090 


146 


4140 


970 


640 


155 


RM^/RM^ 


0.55 


1.18 


0.43 


0.45 


0.35 


0.17 










magnetic turbulence 


d = 2 R c 


= K 


- explosive 


dynamo 














R c 


1060 


461 


461 


922 


1630 


308 


1380 


541 


839 


279 


A B [kpc] 


0.24 


0.068 


0.079 


0.094 


0.19 


0.012 


0.36 


0.31 


0.2 


0.015 


firms b"G] 


4.75 


4.39 


9.44 


10.3 


4.7 


4.06 


5 


3.75 


2.77 


4.09 


RMS [m -2 ] 


443 


155 


762 


1250 


452 


49.3 


659 


247 


131 


54.9 


RM*/RC P S 


3.04 


4.27 


1.57 


2.32 


2.39 


0.51 











Table 1. Application of the model to cool cores of the Hydra A (Hydr.), Centaurus (Cent.), Cygnus A (Cygn.), A1985, A2597, 
3C31, Perseus (Pers.), A85, A2199 and Virgo clusters of galaxies. The first part of the table contains the cluster parameters whose 
values were taken from the literature. References and comments on possible biases of the data are given in Sect. |4] and in detail 
in Appendix The second part describes the expected hydrodynamical turbulence. The third, fourth, and fifth parts describe 
the expected magnetic turbulence for the scenarios of magnetic flux ropes (d = 1, R c = 20), flux sheets (d — 2 R c — 35), and 
flux sheets in the explosive dynamo scenario (d — 2 R c — R*), respectively. Numbers are given here to three digits accuracy, 
irrespective of the fact that the real uncertainties of most quantities are much larger. 



The likely energy sources of the expected turbulence are 
buoyant radio bubbles from the central galaxy, which can inject 
turbulence with the right amount of power, and also on length 
scale, that fit the observed magnetic correlation scales well. 

Motivated by these indications of the physical processes in 
cool core clusters, we developed a steady state scenario for the 
hydrodynamical and magnetic turbulence, assuming that the 
turbulent feedback from the central radio source compensates 
for the cool core radiative losses. This scenario predicts turbu- 
lent length scales and energy densities that can be compared to 
observations. 

The hydrodynamical turbulence could have been probed 
with the Astro-E2 satellite mission which unfortunately failed. 
We hope that successor missions will permit the determi- 
nation of the gas velocity dispersion JSunvaev et alJ l2003t 
llnogamov & Sunvaevl2003l iBrtiggen et all2005al) . At present, 
the magnetic turbulence can be tested with the existing Faraday 
rotation measurements of extended radio sources in cool core 
clusters. The Faraday dispersion is a combined measure of the 
magnetic field strength and correlation length. 



For a sample of prominent cool core clusters, we calcu- 
late the expected hydro- and magnetic-turbulence, and the pre- 
dicted Faraday dispersion. In cases of existing Faraday mea- 
surements, we find that our estimates reproduce the observed 
magnitude of the dispersion over roughly two orders of magni- 
tude in RM (or four in (RM 2 )). On average, our predictions for 
the models based on the flux rope and the fluc tuating dynamo 
JSubramanianll 1 9991 Isubramania n et aljl200fj|) are a factor of 
two higher than the measurements, which is not too surprising, 
since the actual data sample larger cluster radii, whereas our 
estimates are aimed at a radio source at the cluster centre. The 
predictions for the Faraday dispersion in the explosive dynamo 
scenario are a factor of two belo w the observation as shown by 
ISchekochihin & Cowlevl J2006I) . which is more severe, since 
geometrical considerations will only enlarge this discrepancy. 
However, the description of the saturated state of this dynamos 
is currently more a guess than a precise estimate. Thus, also 
this model should be regarded as in agreement with the data. 

Very little fine tuning went into our model, but this is likely 
more a coincidence than a proof of the investigated scenario. 
There are various places in the calculation where factors of or- 
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der one could have been introduced. One example is the as- 
sumption of the bubble radius being exactly the turbulence in- 
jection scale. This assumption will require verification by nu- 
merical simulations of the hydrodynamics, which is beyond the 
scope of this initial investigation. Another uncertain area is the 
adopted parameters of the dynamo theory, which are also not 
fully settled in the literature. 

Nevertheless, we have shown that the straightforward ap- 
plication of current concepts of turbulent magnetic dynamos 
in combination with the emerging picture of cool core stabil- 
isation via heat injection due to the dissipation of turbulence 
seeded by radio galaxy feedback leads to expectations for the 
Faraday rotation signal that match well with the observations. 
This is the case for a variety of galaxy clusters spanning two 
orders of magnitude in their X-ray luminosities. Therefore, our 
picture of cluster cool core heating by radio galaxy feedback 
in combination with ideas about the properties of the saturated 
state of non-helical, small-scale magnetic dynamos passed a 
critical test. Our scenario provides a number of testable predic- 
tions, especially about the level of fluid turbulence in individ- 
ual clusters. We hope that future measurements of the plasma 
properties in galaxy clusters will even permit discrimination 
between the different small-scale dynamo models. 

Acknowledgements. We acknowledge stimulating discussions with K. 
Subramanian and A. Schekochihin on the theory of small-scale dy- 
namos. We are thankful to R. Laing for providing us with the results 
on 3C31 prior to publication, and to V. Springel and E. Churazov for 
comments on the manuscript. 

Appendix A: Individual clusters 

In this Appendix the observational data is discussed. This col- 
lection of literature values f or cool core cluster parameters 
made use of the work done by Pfrommer & EnBlin1 (E()Q4al) . 

A. 1 . Hydra A cluster 

The gas parameters fo r the Hydra A cluster are taken from 
iMohr & Evrardl dl997l) . and the X-ray luminosity is provided 
bv llkebe et alJ dl997lh The first detailed RM map of Hydra A 
was published bv lTavlor & Perlevl(ll993l) . which exhibited very 
large RM values, especially for the south lobe of Hydra A. 
The data was reanalysed using th e PACERMAN algorithm 
( Dola g et aljEool IVogt et alJl20f)i . which removed all areas 
in the northern lobe, and many of the southern lobe with ex- 
tremely high RM values, indicating that those were observa- 
tional artefacts. Since the southern lobe RM data seem to be 
still contaminated by observational artefacts, as statistical tests 
for anti-correlations of gradients of polarisation angles and RM 
values reveal fenBlinetalJ EoCB). we use only the disp ersion 
of the northern lobe, as reported by Vogt et al. (2005). The 
same RM map was analysed bv lVog t & EnBlin (2005) using a 
maximum-likelihood power spectra estimator (see Fig. re- 
vealing central magnetic field strength of B rms as 7 fiG with an 
autocorrelation length of Ag ~ 3 kpc. 

We note that model expectations and observational val- 
ues of Faraday dispersion, magnetic field strength and auto- 
correlation length agree well. Since the Hydra A north lobe is 



known to approach the observer (and the radio emission ex- 
tend to larger radii), the observational RM dispersion should be 
slightly biased to lower values. 

A.2. Centaurus cluster 

The centra l electron density and the cool core radius are 
take n from IMohr et alJ Jl999l) 14 . and the central temperature 
from White (2000) 15 . The observed Faraday dispersion for the 
Centaurus cluster is derived from the embedded radio source 
PKS 1246-410 bv lTavlor et all J2002I) . The published Faraday 
rotation map reveals several isolated patches of extreme RM 
values, and the histogram of RM values exhibits multiple peaks. 
A statistical analysis of the alignment of RM gradients and po- 
larisation angle gradien ts reveals a s ignificant anti-correlation 
between the two (EnB lin et alJl2003l) . which is also an indica- 
tion of observational artefacts in this map. Thus, the reported 
RM dispersion may be too large. However, an order of magni- 
tude agreement between observed and expected RM dispersion 
is found. 

Tayl or et alJ (2002) report a /?M-correlation length of 
isM ~ 0.7 kpc and a field strength of 9 //G, which is signifi- 
cantly different to our expectations. The different length scale i s 
not surprising, since usually Arm > (EnBli n & Vogl 2003). 
The larger observationally derived magnetic field strength may 
be a result of the very small core radius and/ or the possi- 
ble ov erestimate of the RM dispersion reported in lTavlor et alJ 
(2002). 

Thus, the X-ray and radio observational data of the cool 
core of the Centaurus cluster would need further improvements 
before conclusive statements about an agreement or disagree- 
ment of our model and observations can be made in this case. 
Both the observed and the expected RM dispersion could be 
subject to re-adjustments. The complex morphology of this 
cool core probably causes the difficulties singling out a 'real' 
core radius and central electron density. 

A.3. Abell 1958 cluster (3C295) 

The electron density, central tempe rature, cool core ra dius, and 
dispersion measure are taken from Allen et aO (l200ll) and the 
cluster luminosity from I Allen e t al. (2003). 

The agreement of predicted and observed Faraday disper- 
sion is very good. Furthermore, the rather high predicted mag- 
netic field strength of 10 - 20 fiG is also found observationally, 
since lAllen et alJ J200ll) report 14 fiG. 



14 Using the extremely small core radius given in lTavlor et al]l2002T) 
leads to a cool core fractional luminosity o f only 3%, which is obvi- 
ously much smaller than the 40% observed 1 Peres et al. 1998). 

15 Lower centr al temper atures have been reported by 

ISanders & Fabianl <2002l) and iTavlor et alJ J2002l) . but since they 
were not used for the electron density estimates, and since our mag- 
netic field and RM predictions are quite insensitive to the temperature, 
we do not use them. 
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A.4. Cygnus A cluster 

We use the electron density, temperature and cool core radius 
reported in iMohr et alJ 0999I) . We measure the RM disper- 
sion to be 1200 m~ 2 from the map published by Dreheret al. 
(1987) after removing pixels from the map with extraordinary 
large RM jumps, which we associate with observational arte- 
facts (n7r-ambiguities). Similar to the case of Hydra A, we see 
an asymmetry in the RM dispersion between the two lobes, 
which is a consequence of the different depth of the lobes (th e 
Laing-Garrington effect: LaingOHi Harrington et al. 1988). 
Unfortunately, the two radio lobes are located well beyond the 
cool core radius, which is even more true for their polarised re- 
gions used for RM measurements. Since we expect a declining 
magnetic field strength as a function of radius, the reported RM 
dispersion is likely smaller than the one that would be mea- 
sured from a radio source located in the centre of the cool core, 
as discussed in Sect. l3.7l and illustrated in Fig. [2] The predicted 
central RM dispersion is a factor of two smaller than the pe- 
ripherally measured one. With our current understanding of the 
geometry, there is no apparent conflict between these numbers. 



Compton based magnetic field estimates lower by a factor of 
two ( EnBlin et all 19991 lEnfiliil2004l). 

As stated before. [Fabian et alJ J2003"r3) argued for a large- 
scale viscosity in the Perseus cluster cool core of at least 
4 • 10 27 cm 2 /s, which is in good agreement with the 1.6 ■ 
10 28 cm 2 /s predicted in our 1 - d scenario. 

A.6. Abell 2597 

This cool core galaxy clust er has been observed with Chandra 
by McNamar a et alJ (120011). The obse rvation has been reanal- 
ysed recently by Pollac ket alJ (I2005I) . The gas parameters are 
take n from this reanaly sis. The total luminosity used is taken 
fro mlDavid et alJ (Il993|) . 

Ipoflack^UilT O005I) study also the polarisation properties 
of the central radio source in this cluster. They determine a RM 
dispersion of 1080 rad rrT 2 , comparable to our expectations. 
They conclude that the cluster magnetic field has a minimum 
magnetic field strength of 2.1 fiG. Our model predicts a higher 
magnetic field strength of 7-12 uG wh ich is in agreement with 
the lower limit of Pollac k et alJ J2005h . 



A. 5. Perseus cluster 

We take the cent ral electron density, te mperature and the cool 
core radius fro mlChurazov et all J2003I) and the cluster X-ray 
luminosity from D avkTeTal] i 1 993b" 

The Perseus cluster cool core exhibits diffuse radio emis- 
sion, a radio mini-halo. One can apply the classical minimum- 
energy arguments to the radio halo data to search for the 
minimum of the sum of the magnetic and relativistic par- 
ticle energy densities necessary to reproduce the radio syn- 
chrotron errmsi£ri i _arid_assume proton to electron 
ratio k p . Pfrommer& EnBlin (2004b) report a classical mini- 
mum field strength of 7.2^-^fiG (assuming k p = 1), where 
the confidence interval is given by the requirement that the 
energ y density should be within one e-fold from the mini- 
mum. Pfrommer & EnBlin (2004b) also develop and apply the 
hadronic minimum-energy criteria to the Perseus mini-halo. 
This criterion assumes that the relativistic electrons are in- 
jected by hadronic interactions of a relativistic proton popu- 
lation, which usually dominates the relativistic energy budget. 
No proton-to-electron factor k p has to be assumed in this case, 
since the physics of the hadronic interaction determines this ra- 
tio. Applied to the Perseus radio halo, a very similar central 
field strength of 8.8^ 3 4 8 juG was found. Our cool core model 
predicts a field strength of 7 - 13yuG, which is in agreement 
with these findings. A significant lower central field v alue of 
about 1-3 j-iG was reported by ISanders et alJ ( l2005h . based 
on an Inverse Compton interpretation of a hard photon compo- 
nent in the X-ray spectra of the cool core region. If the hard 
photons flux is due to another physical mechanism, the derived 
value will become a lower limit. If the Inverse Compton na- 
ture of the flux could be verified, this will be a very good field 
estimate in the case of a non-intermittent field distribution. If 
there is significant magnetic intermittency, spatially inhomoge- 
neous electron cooling can produce an anti-correlation between 
fields and relativistic electrons, which can easily lead to Inverse 



A.7. 3C31 

The cool X-ray environment of the extended radio source 3C3 1 
is associated with a group of galaxies. X-ray measurements 
suggest that this group shows properties similar to a cool core 
of a galaxy cluster. Therefore, it seems to be a good case to test 
our model in a differen t situation . The p arameters for the gas 
were taken from Hardcastle et al. (2002). However we use the 
total luminosity from Komossa & Bohringer ( 1999). 

Detailed RM measurements have been carried out by Laing 
et al. (in prep.). They determine the RM dispersion to be 25 rad 
m -2 (Laing, private communication). 

A.8. Abell 85 

We use the electron den sity, temperature, c ool core radius, 
and luminosity reported in lMohr et al.l J 1 9991) . and lDavid e t al. 
( 1993J). To our knowledge, no RM measurements or magnetic 
field estimates are available for the centre of this and all fol- 
lowing clusters. Thus, our RM results are predictions. 

A.9. Abell 2199 

The parameters o f the electron density profile as reported in 
IMohr et all dl999l) were used. Th e temperature of the cool core 
are taken fromlVoigt et al.l (12002) and the total luminosity from 
iDavid etalJ(ll993l) . 

A .10. Virgo cluster 

The Virgo clus t er wa s observed with XMM-Newton by 
Matsushita et al ] d2002l) . The gas parameters are taken from 
their a nalysis. However the total luminosity of IDavid etafl 
(1993) is used for our calculations. The large difference be- 
tween the cool core luminosity and the luminosity within the 
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(larger) cooling radius is due to the extremely shallow density 
profile (J3 ~ 0.5) of Virgo. However, this should not affect our 
RM estimates too much, since only a (3 ~ 0.25 should lead to 
problems, as was argued in Sect. 13.71 
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